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We propose a novel diagrammatic method for computing transport coefficients in relativistic quan- 
tum field theory. The self-consistent equation for summing the diagrams with pinch singularities 
has a form of a linearized kinetic equation as usual, but our formalism enables us to incorporate 
higher-order corrections of the coupling systematically for the first time. Furthermore, it is clarified 
that the higher-order corrections are nicely summarized into that of the vertex function, spectral 
function, and collision term. We identify the diagrams up to the next-to-next-leading order correc- 
tions in the weak coupling expansion of (f> 4 theory, which is a difficult task in kinetic approaches and 
other diagrammatic methods. 



I. INTRODUCTION 

The experimental results of heavy- ion collisions at rela- 
tivistic heavy ion collider (RHIC) brought remarkable re- 
sults [1-4]. Among others, elliptic flow in peripheral col- 
lisions suggests that the ratio of the shear viscosity to the 
entropy, 77/s, is smaller than the ever observed in other 
systems [4, 5]. Transport coefficients such as the shear 
77 and the bulk £ viscosities reflect dynamical properties 
slightly away from the thermal equilibrium, while pres- 
sure, entropy, susceptibility, etc., do static properties. In 
the kinetic theory, the shear viscosity is proportional to 
the mean free path of quasiparticles of a long life; a short 
mean free path corresponds to strongly interacting mat- 
ter [6], which is the reason why the quark gluon plasma 
(QGP) at RHIC is called a strongly-coupled quark gluon 
plasma (sQGP) [2]. In any event, the RHIC experiments 
and the subsequent analyses prompted a great interest 
in the relativistic quantum field theory for the transport 
coefficients or non-equilibrium dynamics in general. 

Formally, the transport coefficients are expressed by 
the so called Kubo formula [7] in the framework of linear 
response theory even for the relativistic system; e.g., the 
shear and the bulk viscosities are given by 

r? = lim -Im / d 4 « lu '^(i)([^-(z),^(0)]) , (1) 

C= lim - Im / d 4 xe luJt i9(t)([P(x),V(0)}) , (2) 
>o u> J 

respectively, where V(x) = — T%(x)/3 and Kij{x) — 
Tij{x) + gijV(x). T^ v (x) is the energy-momentum ten- 
sor. We work in Minkowski space with a metric, = 
diag (1,-1,-1,-1). Although concisely expressed by the 
Kubo formula, the computation of the transport coeffi- 
cients is not a simple task in practice. 

Lattice QCD simulation in principle provides us with a 
powerful nonperturbative method, and it is natural that 
there are attempts to apply it to calculate transport co- 
efficients at vanishing chemical potential [8-10]. How- 
ever, one must notice that the calculation of the energy- 
momentum tensor on the lattice is not by far straightfor- 
ward because a more careful analysis would be necessary 



of the low frequency structure of the spectral function 
and the relevant relaxation time scales which are to be 
extracted from the imaginary-time correlation function, 
as warned in [9]. 

Recently, calculations relying on gauge/gravity corre- 
spondence have become popular, in which the retarded 
correlation function can be calculated as the absorbed 
cross section of the black hole [11, 12]. One of the 
most important results from such an approach is that 
the ratio of the shear viscosity to the entropy density in 
N = 4-supersymmetric gauge theories is as small as 1/47T 
with infinite number of colors and strong coupling, and 
this value is conjectured to be a universal lower bound 
[11, 12]. 

Even with the recent development of the fancy meth- 
ods for calculating the transport coefficients, the most 
reliable method with the sound basis in the relativis- 
tic quantum field theory would be diagrammatic ones 
based on the loop expansion. One of the most important 
and difficult part of such a method is, however, how to 
implement a systematic way to deal with the so called 
"pinch singularity" , which would make naive perturba- 
tion theory based on the loop expansion breaks down. In 
fact, this difficulty could be overcome by a resummation 
of specific ladder diagrams. Although restricted to the 
leading order of a coupling constant or large- N expan- 
sion, several methods for such a resummation have been 
proposed [13-17], and their outcome has been shown to 
be equivalent to the results as obtained using the rela- 
tivistic kinetic or Boltzmann equation within the leading 
order of a coupling constant [13, 17]. We here remark 
that direct applications of relativistic kinetic theories are 
also done to compute the transport coefficients in the 
hadronic [18], quark-gluon plasma [19, 20] and the 'semi- 
QGP' phases [21]. 

The merit of the diagrammatic methods based on the 
loop expansion over the direct use of the Boltzmann 
equation is the potential ability to take into account the 
higher-order effects of the coupling constant systemati- 
cally, in principle. As far as we know, there were, how- 
ever, only few works which were able to examine the next 
leading order for a relativistic system, still being confined 
to a scalar theory [22, 23]. The major technical difficulty 
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for going beyond the leading order of the coupling con- 
stant is to classify diagrams order by order in the cou- 
pling constant, because the loop expansion no longer cor- 
responds to the coupling expansion at finite temperature 
and/or density. 

The purpose of this paper is to give a formulation for 
a novel resummation method in which transport coef- 
ficients can be calculated systematically for any order 
of the coupling constant beyond the leading order for 
a relativistic system. Our method is motivated by the 
Eliashberg theory for a nonrelativistic system [24] like 
a Fermi-liquid, although it relies on a cumbersome ana- 
lytic continuation inherent in the imaginary-time formal- 
ism adopted by Eliashberg. We extend and reformulate 
the resummation method developed by Eliashberg to rel- 
ativistic quantum field theory at high temperature. One 
of our ideas is to employ the real-time formalism to avoid 
the complicated analytical continuation, and thus make 
the calculational procedure transparent. The real-time 
formalism is also found convenient to identify diagrams 
corresponding to the pinch singularity, as will be dis- 
cussed in Sec. IV. 

The paper is organized as follows: In Sec. II, we briefly 
review why the difficulty arises when computing the 
transport coefficients. In Sec. Ill, we review the real-time 
formalism in the R/A basis that is useful to identify the 
dominant diagrams in the Green function corresponding 
to the transport coefficient. In Sec. IV, the computation 
method for transport coefficients are discussed, and we 
show that the resummation equation has a form similar 
to a linearized Boltzmann equation. In Sec. V, we will 
conclude the paper with some outlook. 



II. PINCH SINGULARITY AND 
QUASIPARTICLES 

The difficulty in computing transport coefficients in 
Eqs. (1) and (2) arises from the infrared limit, cj — > 0, 
of the Green function corresponding to a long time scale. 
In such long time scales, a lot of microscopic scatterings 
occur. In a diagrammatic method, such scatterings are 
expressed by multiple loop diagrams. The question is 
what kind of diagrams are dominating and should be re- 
summed. A particle excitation of a long life is called a 
quasiparticle, and as long as the quasiparticle descrip- 
tion is valid for the system under consideration, thermal 
excitations of the system will be well described solely 
by the quasiparticle excitations. Such an excitation is 
expressed by a product of the retarded and advanced 
Green functions of the quasiparticles, which product is 
found to become divergent in the infinite lifetime limit. 
This divergence or singularity is nothing but the pinch 
singularity that we mentioned in Introduction. 

To see this more concretely, let us take the </> 4 theory, in 
which the one-particle retarded and advanced propagator 



read 



k 2 - m 2 — Re II R (k) — i Im II R (k) (3) 
= D A *(k), 

where to is the mass of the scalar particle, and II R (k) 
is the retarded self-energy. These propagator can be ap- 
proximated by a sum of poles corresponding to quasipar- 
ticle excitations and thus we have 

D R (k) ~ V ~ z n{k) . . 

DA ^^ k °-e n (k) k \ ln (k)> (5) 

where e n (fe) is the energy of the quasiparticle, z n (k) the 
renormalization function, and 7«(fe) > the damping 
rate. At weak coupling and high temperature (to <C AT), 
e n (fe) ~ ±yjk 2 + rrij,, with the thermal mass defined in 
Eq. (89) bellow, and z n (k) ~ l/(2|e„(fe)|). The damp- 
ing rate owing to 2 — > 2 scattering is of order A 2 T. In 
the naive perturbation theory, one would encounter a 
product of the retarded and the advanced propagators, 
D R (p + k) D A (k) , which has the following anomalous be- 
havior at small p under the quasiparticle approximation: 

D R (p + k)D A (k) 

where v n = de n (k) /dk is the velocity of the quasiparticle. 
Equation (6) shows that when the damping rate vanishes, 
i.e., 7n(fc) = 0, the product diverges at p = 0; this is the 
pinch singularity as mentioned above. And one sees that 
the product of the retarded and the advanced propagator 
with the common momentum must be summed over to 
have a sensible result. 

The need for resummation can be seen even apart from 
the quasiparticle approximation. The product of the 
propagators at p = exactly becomes 

lim D R (p + k)D A {k) = / ( *) , (7) 
p^o v -2Imi7 fi (fc) 

where p{k) is the spectral function defined by 

p(k) = 21m D R (k) 

= -2ImJ7 fl (fc) (8) 

(fc 2 - to 2 - ReTJfl(fc)) 2 + (lmiI R (fc)) 2 ' 

One sees that Eq. (7) diverges when the imaginary part 
of the self-energy vanishes. Therefore, we have to employ 
the dressed propagator, Eq. (3), to avoid the singularity. 

III. REAL-TIME FORMALISM 

In this section, we briefly review the real-time formal- 
ism [25]. For simplicity, we consider the </> 4 theory. The 
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FIG. 1. A complex-time path in the real-time formalism. 



Lagrangian has the form, 



X 

4!' 



(9) 



The real-time formalism is formulated on a complex-time 
path shown in Fig. 1. The fields on C\ and C2 are called 
type-1 field, <j> 1 , and type-2 field, </> 2 , respectively. The 
partition function is 



Z12 = / V^Vcf 1 exp 



1 



d 4 xd 4 



-i f dS(V(^)-V(0 2 )) 



(10) 



Although there are contributions from the paths C3 and 
C4 in general, they are factorized in the limit of t j, tf — > 
00, and are irrelevant as long as correlation functions are 
concerned [25] . The free propagators in momentum space 
read 



D$ = (1 



D 



D 12 _ 



n{k°))D F (k) + n(k )D* F (k) , 
e° k ° (n(k°) + 9(~k )) (D F (k) + D* F (k)) , 
e-° k ° (n(A°) + 6(k )) (D F (k) + Dp(k)) , 
Df = (1 + n(k a ))D F (k) + n(k°)D F (k) , 

where the a is the parameter characterizing the time path 
in Fig. 1, D F (k) is the free Feynman propagator defined 
by 



D F (k) = 



(12) 



For later convenience, let us rewrite Eq. (11) in terms 
of the retarded and advanced propagators, 

D F \k) = (1 + f(k )) (-iD*(k)) + f(k°)(iD F (k)) , 

Df{k) = e^°f(k )(-zD§(k)+iD F (k)) , 

Df{k) = e-° k ° (1 + /(fc )) (-iD^(fc) + iD F (k)) , 

Df(k) = f(k°)(-iD*(k)) + (1 + f(k°))(iD F (k)) , 

(14) 

where the free retarded and advanced propagators are 
defined by 



-1 



fc 2 



ek° 



D F (k) 
and, f(k°) is defined as 



-1 



k 2 — m 2 — iekP 



1 



(15) 
(16) 

(17) 



eP k ' - 1 

These propagators satisfy Dp(—k) = Dp* (k) = Dp(k), 
which relation also holds for the full propagators (see 
Appendix A.) The relation between the Bose-Einstein 
distribution function and f(k°) is n(k°) = f(\k°\). 



A. R/A basis 

In this subsection, we introduce a useful basis called 
U R/A basis" for computing transport coefficients [26, 27]. 
In the original representation introduced by Aurenche 
and Becherrawy [26], the propagators are diagonal, and 
the diagonal components are the retarded and advanced 
propagators. In the present work, we rather employ the 
representation in which the diagonal components are zero 
while the off-diagonal components are not [27] , instead of 
the original R/A basis (we also refer this as R/A basis in 
this paper.) We find that this basis is useful not only for a 
comparison of the real-time correlation function obtained 
in the real-time formalism with that in the imaginary- 
time but also for identification of pinch singularities. 

Let us define a rotation matrix for converting the basis 
Eq. (11) (we refer this basis as "standard basis") to the 
R/A basis in momentum space as 



(b a (k) ^U a a {k)<j> a {k), 



(18) 



and n{k ) is the Bose-Einstein distribution function, 

< k0 ) = ^wv^~v (13) 

with inverse temperature /3 = 1/T. We note that any 
physical observable is independent of the choice of the 
path, i.e., a. 



where the Greek index denotes R or A, and the Latin 
index denotes 1 or 2. We use Einstein notation, i.e., if an 
index appears twice in a single term, once as a superscript 
and once as subscript, a summation is assumed over all of 
its possible values. We also define the metric for the stan- 
dard basis as g a b = g ab = diag(l,— 1); then, the metric 
in the R/A basis is given by g aj 3 = g a bU a a {k)U b p(—k) , 
of which explicit form is given below. The inverse matrix 
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U b a (k) is defined by 5 b = U b a {k)U a a {k). An n-point 
function transforms as a tensor: 

r a /3...(ki,k 2 , ■ ■■) 

= r ab ...{k ll k 2 ,---)u a a {k 1 )u b S) {k 2 )--- , (19) 

£>^-(fci,fe 2 ,---) 

= D ab -(k 1 ,k 2 , ■ ■ • )U a a (k 1 )U p b (k 2 ) • • • . (20) 

In particular, the propagator or two-point function trans- 
forms as 



D ab {k) = U a a {k)U b {-k)D al3 {k) . 



(21) 



In the R/A basis, D a @(k) is chosen so that the diagonal 
part vanishes: 



n^(k\ ( DRR W D RA {k) 

U (k) ~ \D AR (k) D AA {k) 

-iD R (k) 
iD A (k) 



(22) 



The most general form of U a a (k) satisfying Eq. (22) is 
obtained through comparison with Eq. (14): 



TT« (k) - ( U± ^ k ) U1a ^ 



/(fc°)e CT ' c0 c(fc) e^-^" /c(-k) 
f(k°)c(k) l/c(-k) 



(23) 



where c(k) is an arbitrary function. Here, we simply fix 
c{k) to be constant, c(k) = 1, while the a is kept to be a 
free parameter. Then, 



a[ ' ~ \ f(k°) 1 



(24) 



and the inverse is 

rra (k] _ (U R l{k) U R 2 (k) 

u aW ~ \U\ik) U A 2 (k) 



e {P-cr)k a 

(l + f(k ))e-° k ° (l + /(fc°)) y 



(25) 



The metric in the R/A basis is 

g aP = ga b U a a (k)u b (j {-k) = (J J) = 9 al3 



(26) 



which is independent of a. The propagator can be rewrit- 
ten by using the metric g a ^: 



D a0 {k) = ~ig af3 D a {k) 



(27) 



with no sum over a. Here we define the Feynman rule 
for the propagators with arrows: 



-iD R (k), (28) 
-W A (k) = -iD R (-k). (29) 



The propagator is identified with D R (k) or D A (k), if 
the arrow is parallel or anti-parallel to the momentum k, 
respectively. The four-point vertex transforms as 



X a pys(kl, &2) h, foi) 

= X abcd U a a (k 1 )U b ^k2)U%(k 3 )U d s (h) , 



(30) 



where X a bcd = +X i£a=b = c = d=l, and \ a bcd = — X 
if a = b = c = d = 2, otherwise X a b c d — 0. The vertex 
satisfies the energy-momentum conservation, ki + k 2 + 
&3 + fc 4 = 0. After simple calculations, one finds 

X a p~fs(ki,k 2 , fc 3 , ki) 

= A[(l + /(fc?))^(l + /(fc°))^ 

x(l + /(fc 3 °))^(l + /(fc 4 ))^ 

" (/(fc 1 ))^(/(fc 2 °)) 5 ' 3i? (/(fc3 ))^(/(fc4 ))^" • 

(31) 

The first and second terms in RHS corresponds to the loss 
and the gain terms, respectively. The four-point vertex 
generally contains sixteen combinations; however, thanks 
to the symmetry of the </> 4 theory, they reduce to the 
following five vertices: 




= -iXRRRR(ki,k 2 ,k 3 , fc 4 ) = , 



= -i^RRRA(ki,k 2 , k 3 , fc 4 ) 



(32) 



X 

X 
X 



-iA [(l + /(*?)) (l + /(*§)) (l + /(*§)) 

- /(*?)/(*£)/(*:§)' , 

(33) 

: -iXR R AA(ki,k 2 , k 3 , fc 4 ) 



-iX 



(l + /(*?)) (l + /(fc 2 )) 
-/(*?)/(*§)" , 



(34) 



= -iXRAAA(ki,k 2 ,k 3 ,ki) = -iX, (35) 



-iXAAAA(kl,k 2 ,k 3 ,k4,) = 0. (36) 
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At the tree level, vertices with all the same indices van- 
ish, i.e, Xrrrr = ^aaaa = 0. This relation is satisfied 
in full n-point vertex functions; see Appendix A for the 
derivation. We should note, however, that this vanish- 
ment will be no longer maintained for the system out of 
equilibrium, where Kubo-Martin-Schwinger condition is 
not satisfied [28]. 



IV. COMPUTATION METHOD FOR 
TRANSPORT COEFFICIENTS 

In this section, we discuss the computational procedure 
for transport coefficients in the R/A basis, and derive a 
self-consistent equation, which turns out to be an exten- 
sion of the linearized Boltzmann equation. We shall see 
that the R/A basis can naturally lead to a decomposition 
of diagrams into the dominant diagrams including pinch 
poles and others. 



A. One-loop analysis 

We shall start with the one-loop diagram of a re- 
tarded function. Here, we employ the dressed propa- 
gators, D R (k) and D A (k), which do not have a pinch 
singularity. We assume that the relevant operator 0(x) 
to the transport coefficient is a quadratic one. The re- 
tarded Green function of it reads 

G%(p) = J d^xe^emOix), O(0)}) . (37) 

At one-loop level, 



2 j ^ r A% ai (-P>- k >P + k ) 

x D a ^(p + k)D^{-k) ( 38 ) 



where the overall of 1/2 is a symmetric factor, 
and the Fcynman diagram for the vertex function, 
^a3-y(ki; k 2 , k 3 ) corresponding to 0(x) is given by 




/.-, = r ARR( k ^ k 2,k 3 ) 



= rW(h,k 2 ,k 3 ){i + f(k° 2 ) + f(k° 3 )) , 

(39) 



(40) 



and r^ RR (ki,k 2 ,k 3 ) = r^l A (ki,k 2 ,k 3 ) = 0. For ex- 
ample, rW(k u k 2l k 3 ) = k 2 k\ + k\k{ - 2S ij k 2 ■ fc 3 /3, for 
the shear viscosity, where 0(x) cx 7Ty(a;). There are, 
in general, four diagrams shown in Fig. 2 contributing 
to the retarded Green function. The diagram (4) van- 
ishes because the vertex functions are r RRR (ki, k 2 ,k 3 ) = 
r i A ) AA (k 1 ,k 2 ,k 3 )=Q. The diagram (1) is expressed as 



d A k 



(1 + f(p°) + f(k )) 



2 V J (2tt) 4 
x r (0) (-p, -k,p + k)D R (k)D R (p + k) 

x r (0) {p,- P - k,k) 
2 / ^-J( k °)r {0) (-P,- k ,P + k ) 

x r<°> (p, -p - k, k)D R (k)D R (p + k) , 

(41) 



where we have dropped the term proportional to l+/(p°) 
that does not contain poles in the complex A; plane, since 
the vertex and the retarded propagator have no poles in 
the upper complex k° plane; the k integral becomes zero. 
For a transport coefficient, we need the soft-momentum 
limit as 



<t« = lim lim \lmg^ (1 \p) 



pO-i-Op-^O p° 

i r d A k d 

4 J (2tt) 4 dk° 



f( k °) 



(42) 



x (rW(0,fc,-fc)) 2 Re(D R (k)) 2 . 

For the contribution Q R (p) from the diagram (2), one 

finds — The contribution Q R (p) from the 

diagram (3) includes a product of the retarded and ad- 
vanced propagators: 



<r '(P)=2H) 



d A k 



{ , 7 _ )A (l + f(-k )+f(p° + k )) 

x D A (k)D R (p + k)r w (-p,-k,p + k) 

x r {a) ( Pl -p-k,k) 
i f d"k 



if(p° + k°)~f(k )) 



(2tt)* 

x D A (k)D R (p + fc)r (0) (-p, -k,p + k) 
x r^(p,-p-k,k). 



(43) 



Then, the transport coefficient due to Q R (p) be- 



comes 



r(3) 



lim lim —ImQ R ^ 3 \p) 



p°->0 p->o p 



1 



f(k°) 



d A k d 
(2tt) 4 



(44) 



G 




FIG. 2. One-loop diagram for a retarded Green function. The red line corresponds to the product of the retarded and advanced 
propagators. 



If one uses the quasiparticle approximation, p(k) = 
27re(k )S((k ) 2 - El) and lk = -ImII R {E k ,k)/(2E k ), 
<j( 3 ) becomes 



r(3) 



d 3 k 

(27rf 



f{E k ){\ + f{E k )) 



i /r(°)(o,fc,-fc) 

Wk \ 2E k 



(45) 



while and vanish in this limit. At weak coupling, 
the decay width is proportional to A 2 , so that cr'- 3 ' ~ 
1/A 2 . This gives the correct coupling dependence, but 
not the coefficient. The higher-order diagram contributes 
to the transport coefficient in the same order as the one- 
loop diagram. 

Since the Bosc-Einstein distribution function diverges 
at fc° = 0, each (p) includes divergence, although 

the sum of G R (p)'s is finite. In order to avoid such an 
artificial divergence, we introduce an infrared (IR) cutoff 
k. Bosons with very soft momentum can no longer be 
regarded as a particle, so we have to treat it rather as 
wave. This IR cutoff n separates the scale between the 
particle and the wave. Thus, we take the IR scale to be 
k -C m. If one applies the quasiparticle approximation, 
the Ilft(k) become independent of k, because k° > m^> 
k at the on-shell of the quasiparticle. We note that a 
fermion has no such a divergence because of the Fermi- 
Dirac distribution function. 

The merit of the R/A basis is that it enables us to iden- 

tify the diagrams with large contribution to Q R (p). Al- 
though such an identification and decomposition are not 
obvious in the imaginary-time formalism before the ana- 
lytical continuation, the decomposition plays important 
role for the full analysis of the retarded Green function, 
as is discussed in the next subsection. 



Relativistic Eliashberg Decomposition in 
Real-time Formalism 



In the previous section, we showed within the one-loop 
level that the product of the retarded and advanced prop- 
agators with the common momentum gives a large con- 
tribution to transport coefficients. Such a product of the 
propagators also appears in higher-order diagrams, which 
must be also summed over. Eliashberg method [24] is one 



o<i> 



FIG. 3. A retarded Green function corresponding to a trans- 
port coefficient. The dark (green) square denotes the con- 
nected four-point function. The Lines without arrows denote 
lines that can be the retarded or the advanced arrow lines. 



of such a resummation method for non-relativistic sys- 
tems and has been applied to Fermi liquid at low temper- 
ature in the imaginary-time formalism. In this method, 
the four-point function is first analytically continued to 
the real-time domain. Each external leg can be the re- 
tarded or advanced legs, so the four-point function has 
a matrix form. After the analytic continuation, the di- 
agrams for the four-point function is decomposed into a 
set that connects with pinching diagrams and a set that 
does not connects with them. The pinching diagrams 
are summed over by a self-consistent equation, which is 
found to correspond to a kinetic equation. In this sec- 
tion, we extend Eliashberg method to a relativistic case 
using the real-time formalism in the R/A basis, and thus 
derive the corresponding relativistic kinetic equation. 

The retarded Green function shown diagrammatically 
in Fig. 3 is expressed in the momentum space as 



Sg(p) = J dW*-*i0(t)<[0(s),0(O)]) 

_i r d A k r d 4 fc' 

" 2 J J2t^J (2<r 

+ H)T aMa3 (-p - k, k, -k',p+ k')D a ^(p + k')D^(-k')] (p, -p - k', k') , 



i r Ap lQl ("P. ~ k >P + k ) Daia2 (P + k ) DPl02 (- fc ) [(2^) V 4) (k - k')5 a . 2 a -8 h ^ (46) 



where Ta 2 f3 2 p 3 a 3 {—p ~ k, k, —k',p + k') denotes the con- nected four-point function. The first diagram in RHS of 
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Fig. 3 is just the one discussed in the previous subsec- 
tion. Owing to Eq. (27), the product of two propagators 
appearing in Eq. (46) is written as 



D aia2 (p + k)D^^(-k) 

= H) V 1 " 2 g M2 D ai (p + k)D h (-Jfc) 
= -g aia2 g plp2 G aiPl ( P ,k). 



(47) 



= 4 0) (p,*) 

= T (0) (p, -p-k,k) 

x (/(p°) - /(p° + fc )) , 

= 4 0)3 (p,fc)-r(°)(p,-p-fc,fc), 



(51) 



Thus we have 



x G Qlft (p, Jfc) [(2^) 4 <5 (4) (fc - fcW% a 
+ *T ttlft /32Q2 (-p-ft,fc,-ft',p + fc') 



x G a2f:)2 ( P ,k') 



r. 



(0) 



(p, -p - ft', ft') . 

(48) 




(52) 



+ fc = r<™ (p, fc) = r(°) (-p, -fc, p + k) , 



(53) 



= ri° 2 ) (p,fc) = r(°)(-p,-fc,p + fc), 



(54) 



The indices appear in pairs of a and (3 in Eq. (48), so it 
is useful to introduce an index I as the pair of R and/or 
A, i.e., I = (a,/?) = (RA), (AR), (RR), and (AA), then 



d 4 /c 



d 4 fc' 



(2tt) 4 J (2tt) 4 
(27r) 4 (5( 4 )(fc-fc')^ m 
+ iT i m (p,fc,fc')G m (p,A : ') 



4 0) ""(p.o. 

(49) 



where the three-point functions are 




4 0) V*o 



r(°)(p,-p-fc,fc)(i + /(p°) + /(fc )), 

(50) 



= r(°) (-p, -fc, P + k) (f(p° + k a )-f(k )) 



and r^' (p, fc) = r^ 1 (p, k) = 0. At the tree-level, the 
three-point function satisfies 

4%, fc) = (/(p° + k°) f(k°))r^(- P , -k) . (55) 

This relation is generalized to the full vertex function as 

r A3 (p,k) = (/(p° + fc°) - f(k°))n 3 (-p,-k) , (56) 

the proof to which is given in Appendix A. 

As discussed in the previous section, G 3 (p, k) contains 
pinch singularities at weak coupling limit. Therefore we 
treat G 3 (p, fc) separately from other G l (p,k) with i ^ 3 

in the four-point function. Suppose that T^ m (p,k,k') 
is the four-point function that does not include a pair of 
lines of the type G l (p, k"). Then, the four-point function 
obeys the following equation: 



,(0) 



T l m (p,k,k , )=T {A) m(p,k,k') 

.„ ; a ft 
+ (-*) 



rl;)4 T^ n {p^k") (57) 
xG"(p, k")T n m {p,k",k'). 

In order to pick G 3 (p,k) up, we also define 

T*- 5 ' m (p,k,k') as a four-point function that does not 
include G 3 (p 4- k",k"). The full four-point function 




FIG. 4. Self-consistent equation for the four-point func- 
tions: T 3 3(p,k,k'), T %{jp,k,k'), and T^P, k, k'), respec- 
tively. The blue square denotes T x j(p,k,k') and the green 
square denotes T^ s '*j(p, k, k'). 



T s(p, k, k') obeys the following self-consistent equation: 

TMp,k ) k') = ^ 3 3 (p,k,k') 

d k r>\ 3 



T (B) 3 {p,k,k") (58) 



(2tt)4 

xG 3 (p,fc")T 3 3(p,fc",fc') 



Other full four-point functions satisfy 

T i j ( P ,k,k') = T^ i j ( P> k,k') 

+ H) 3 / 0^T^\( P ,k,k") 

/rl 4 h" r rl 4 lc'" 

xG 3 (p,k")T 3 3 (p,k",k'") 

xG 3 ( P ,nr (B)3 j {!> ) r ) fc'), 

(59) 

r 3 ( P ,k,k') = T( B) \(p,k,k f ) 



+ j ~T^\{p,kX) ( 60 ) 
xG 3 (p,k")TM P ,k",k'), 



T 3 j {p,k,k') = T^ jfak,*) 



+ H) 3 / ^r 3 3b,fc,fc") 
xG 3 (p,nr (B)3 >,fc",£0- 



(61) 



The diagrams corresponding to Eqs. (58) to (61) are 
shown in Fig. 4. Inserting Eqs. (58) to (61) into Eq. (49), 



FIG. 5. Vertex corrections to F(p,k). The green square cor- 



responds to T' s ' i(p,k,k'). 

we arrive at a simple form: 
ft-\ % f d k f d k 

Qr{p) 



r*(-p,-k) 



(2tt) 4 J (2tt) 4 
x (/(fc°)-/(P° + fc°))G 3 (p,fc) 

(2^) 4 <SW(fc - k') + iT 3 3 ( P , k, k')G 3 (p, k') 

xr{ P ,k') + K{p), 

(62) 



where 



K{p) = 



2 

-Y 

2 ^ 



d*k 



x 



(2tt)< 
[(27r) 4 ^W(fc-fc')<5 l J 
+ iT {B)l J (p,k,k')Gi(p 7 k') 



r 



(o)J 



(63) 



which does not include G 3 (p, k), so that JC(p) is not en- 
hanced at small p. We have also introduced the effective- 
vertex function shown diagrammatically in Fig. 5, which 
is expressed as 



r( P ,k) = r^(p,k) 



+ H) 3 E/ 



d A k' 



T^\( P ,k,k>) (64 ) 

(0)* 



(2tt)* 



This effective- vertex function includes vertex corrections 
which can be evaluated by loop expansions as long as 
further infrared singularities do not appear, and contains 
both quantum and medium effects; these corrections are 
in the same order in coupling at high temperature. We 
emphasize that this vertex correction is not included in 
the usual kinetic theory, and hence the vertex correction 
derived here describes an effect beyond the kinetic theory. 

Now the transport coefficient a is given by the retarded 
Green function at the p — > limit, 

a = lim lim — lm Q R (p) = a B + <?k ■ (65) 
p°->op->o p v 

Note that these limits, p Q — > and p — > 0, do not ex- 
change in general. We have decomposed the transport 
coefficient into two parts, ctb + &k.i the explicit forms of 
which will be given shortly. 



Let us discuss the two terms one by one. The first term 
is expressed as 



1 



ctb 



d A k 



d A k' 



(2tt) 4 J (2tt) 4 

x P( k ) 
-2Im7T fl (fc) 



r*(-k) 



d 



lmT 3 3 (k,k') 



<9fc° 

{2n)H {i) {k - k') 
p(k>) 



f(k°) 



(66) 



-2Im77 fl (fc') 
where we have used the relation 

p[k) 



r(fe'). 



G 3 (0,fc) = D R {k)D A (k) 



-21m n R (k) 



(67) 



in accordance with Eq. (7). We have also introduced 
abbreviated notations: 



T 3 3 (k,k') = T 3 3 (0,k,k / ), 

r(k') = r(o,k'), 
r*(k) = r*(o,k). 

The second term of RHS in Eq. (66) is defined by 



1 



Ufr = lim lim — Im/C(p) . 



(68) 
(69) 
(70) 



(71) 



As will be clarified in the next subsection, the physics 
content of the <jb can be nicely given in terms of a Boltz- 
mann equation. At weak coupling, as becomes much 
larger than ujq because of the pinch singularity. In the 
following, we focus on as- 

Let us introduce the full vertex function, -T(fc), by 



r(k) = F(k) 



d A k' 

W) 1 



lmT 3 3 (k,k r )- 



P^) r(k') 
-2lmn R (k>) y n 

(72) 



which is found to obey the following self-consistent equa- 
tion on account of Eq. (58): 



r(k) = r(k) 



d 4 k' 



-ImT (S)3 3 (fc,fc')- 



p(k') 



-21m II R (k') 



(73) 



The diagrammatical representation of this integral equa- 
tion is given in Fig. 6. 

Now ' 3 (p,k,k') consists of two parts: One is 

T^'' ' 3 (p,k,k') that does not include G l (p,k), and the 

other G l (p,k). The latter part contains 3 (p, k, k'), 
which is proportional to f(p° + k'°) — f(k'°), so that 

it vanishes as p — > 0. Therefore, T {B) 3 {k,k') = 

3 (k,k'). The dominant transport coefficient as is 



> > J> 



FIG. 6. The self-consistent equation for the vertex function. 
The blue and green blob corresponds to the full and renor- 
malized vertex function. The green square corresponds to the 
collision term. 



now expressed as 



OB 



1 

"2T 



d 4 k 



p(k) 



(2tt) 4 -2Im7T fl (fc) 

f(k ){i + f(k°))r*(k)r(k), 



(74) 



where we have used the relation r*(k) = r*(—k). 

Equations (64), (73), and (74) are our main result. 
What we have done is rewriting diagrams in a useful and 
physical form. This is just a rewriting of diagrams, so 
that there is no approximation in the sense of the dia- 
grammatic expansion (although nonperturbative contri- 
butions such as instantons are not included.) The n R (k), 

T {A) 3 (k,k'), and f(k) can be expanded by loops as 
far as further infrared singularities appear. The spectral 
function is obtained through II R (k). One has to choose 

the loop diagrams for n R (k) and 3 (k, k') so as to 

satisfy the symmetries of the action, i.e., Ward-Takahashi 
identity, as is discussed in Appendix C 2. 

Let us estimate the order of <jb at weak coupling. The 
leading contribution comes from the peak of the spectral 
function. The residue of the pole corresponding to the 
peak is of order 1, and Im7T/j(/c) is of order A 2 at the 
pole, so the first part in Eq. (74) gives a contribution of 

order 1/A 2 . The four-point function, ImT^' 3 (k,k r ), in 
Eq. (73) is related to the squared scattering-amplitude, 
which is of order A 2 from 2 — > 2 scattering. This is the 
same order as the imaginary part of the self-energy, so 
the full vertex -T(fc) is the same order as r(k). As we will 

see in the next subsection, ImT'" 4 ' 3 (k, k') and Im n R (k) 
are related to the collision term of a Boltzmann equation, 
which has the same order of the coupling constant. As 
a result, the is estimated as of order 1/A 2 at weak 
coupling. This is consistent with the result in kinetic 
theory where the shear viscosity is proportional to the 
inverse of the transport cross-section of order A 2 . We 
note that the <jb diverges at the zero cutoff limit, k = 0, 
for bosons, because f(k) has the pole at k° — 0, while 
fermions do not. The divergence will be cancelled by 
adding ok,. 
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C. Linearized Boltzmann Equation 

Here we discuss the relation between our diagrammatic 
method and the linearized Boltzmann equation. We will 
show that Eq. (73) has the form of a linearized Boltz- 
mann equation. In the leading order of the coupling con- 
stant, it is known that the linear equation, Eq. (73), with 
the quasiparticle approximation is reduced to a linearized 
Boltzmann equation in scalar theory [13] and QED [17]. 
Beyond the leading order, nothing has been known about 
the relation between diagrammatic method for Kubo for- 
mula and kinetic equation. Although the Boltzmann 
equation is the equation for the on-shell (quasi-) par- 
ticles, the spectral function includes not only the quasi- 
particle peak but also multi-particles state spectrum. If 
one wants to derive the Boltzmann-like equation, which 
contains the scattering process of the quasiparticles, one 
needs to decompose the spectrum function into the quasi- 
particle part with a distinct peak and others, and rewrite 
the equation in terms of the quasiparticles. In this paper, 
we shall not do such a decomposition; nevertheless, we 
shall show the linear equation has the similar property 
to that of linearized Boltzmann equation. To derive a 
linearized collision operator in our formalism, we define 
ip(k) by 



<p(k) 



r(k) 



-2Imi7 fl (fc) 



(75) 



following [13]. Then, Eq. (73) is cast into the following 
form 



—r(k) = cm 

where C is defined by 

dp(k) = -2 ^Im n R (k)ip(k) 



(76) 



i r d A k* 



2k° J (2tt)' 



p(k')lmT {A)3 3 (k,k')if(k'). 



(77) 



We note that ip(k) is an odd (even) function under k <H> 
— k, if -T(fc) is an even (odd) function, since Im II R (k) is 
an odd function. Equation (76) has the same form as 
the linearized Boltzmann equation [29]. In fact, C can 
be identified with the linearized collision operator: The 
equivalence between Eq. (76) and the linearized Boltz- 
mann equation in the leading order of the coupling con- 
stant is shown in Appendix B. 

Here we introduce an inner product as 



(cpi,cp 2 ) = 



d A k 



W(k)vUk)<P2(k) 



(78) 



where W{k) = fc°(l + f (k )) f (k°) p(k) > is the weight 
function, and (pi and if2 are arbitrary functions of k^. 
We note that W(k) at k — is finite, although f(k°) 



has a singularity at fco — 0. Using this inner product, we 
can rewrite the transport coefficient, Eq. (74) as 



* B = -(S,£- 1 S) 



(79) 



where S = r/(2k°). Note that S must be orthogonal to 
the eigenvectors of C with zero eigenvalues corresponding 
to conserved charges; otherwise £ _1 is not well-defined. 
The inverse of the collision operator can be expressed 



die 



-tc 



(80) 



where we have assumed that the eigenvalues of C are 
nonnegative. This is necessary for the stability of the 
system. Using Eq. (80), Eq. (79) is rewritten to a Kubo 
formula for the quasiparticle, 



ob = 



dt(S(0),S(t)) , 



(81) 



where S(t) — exp(— 1£)6>(0). We rewrite a Kubo for- 
mula in field theory to that for relevant quasiparticles. 
Equation (81), a semiclassical formula, i.e., no field op- 
erator is present in Eq. (81). The collision operator and 
the effective vertex are calculated in thermal field theory. 
We note that the thermal weight W(k) does not contain 
f{k°) but /(/c°)(l + f(k )), because the thermal average 
is taken for a two-point correlation function. If the S is 
an eigenstate of C with an eigenvalue, 1/t, 



T 



(5(0), 5(0)) 



(82) 



The inverse of the eigenvalue r can be identified as a 
relaxation time. 

Let us estimate the shear viscosity using Eq. (82) in the 
quasiparticle limit, where, the thermal weight function 
becomes 

W(k) ~ 2n[S(k° - \k\) + 5(k° + |fe|)]/(|fe|)(l + /(|fe|)) , 

. (83) 

and the vertex function for the shear viscosity tt ij in the 
leading order is 



k° 



k l k 3 -5 13 - 



Then, the shear viscosity is evaluated as 
WT { ,nij) 5 ' 



V 



(84) 



(85) 



where s = 27r 2 T 3 /45 is the entropy density in the free 
theory. Using P = 4s for free massless particles in the 
leading order, we obtain the relaxation time as 



5 77 
4P 



(86) 



This is 25% larger than that obtained in Ref. [30], al- 
though Eq. (85) is a parametrically similar to that given 
in Ref. [30] . The difference comes from the semi-classical 
calculation in our formalism and the field theoretical cal- 
culation in Ref. [301. 
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FIG. 7. The diagram contributing to T (A) s(k,k') in the 
leading order, which corresponds to 2 — > 2 scattering. 



D. Leading and higher orders 

In this subsection we show how the contributions of 
the higher as well as leading-order terms to the transport 
coefficient appear in our formalism. The transport coeffi- 
cient (7b in Eq. (79) consists of the vertex function S, the 
collision term C and the inner product that contains the 
spectral function p(k). The contribution of the leading 
and higher order terms are nicely summarized into these 
terms. The merit of our formalism is that the contribu- 
tions of each term can be systematically estimated using 
loop expansions. To show the details of this statement, 
let us consider the shear viscosity as an example. We 
will find the shear viscosity is expanded by the coupling 
constant A as 



^ = 1 + CrN/A + (c' 2 In 1/A + c 2 )A + 
Vo " ~ 



(87) 



Here the numerical calculation in the leading-order gives 
t] ~ 3033.54T 3 /A 2 [13, 22]. 

First, let us begin with clarifying how the A~ 2 - 
dependence of 770 is obtained in the leading order within 
our formalism, although a biref discussion on this matter 
was given in the end of subsection IV B. The diagram of 

T( A ) 3 (fc, k') in the leading order is shown in Fig. 7 that 
corresponds to the 2 — > 2 scattering. Then the squared 
scattering amplitude is proportional to A 2 , and thus the 

collision term is estimated as £ ~ 3(k,k') ~ A 2 . 

In this calculation, the quasiparticle approximation with 
a vanishing mass is employed in the spectral function, 
p(k) = (27r)sgn(fc°)(5(fc 2 ), because hard momenta of or- 
der k ~ T contribute to the leading order, where the 
mass can be neglected. Formally, the diagram in LHS 
of Fig. 7 also contains 1 — >• 3 decay and 3—^1 fusion 
processes; however, they are found to be higher orders of 
the coupling because of the phase space restricted by the 
enegy-momentum conservation. The contribution of S is 
of order one. As a result, we have r) ~ S 2 /C ~ 1/A 2 , 
which coincides with the result given by the Boltzmann 
equation. In fact, the equivalence between the diagram- 
matic method in the leading order and the linearized 
Boltzmann equation generally holds, as is proved in ap- 
pendix B. 

Next, let us estimate the contribution of the next lead- 
ing order (NLO) to the shear viscosity in Eq. (87), which 
contributions are found to come from the correction to 
p(k), while the loop corrections to S and C, at least of 



D*£x>Da>U> 



FIG. 8. Vertex corrections up to of order A 2 



order A, are higher-order corrections. The spectral func- 
tion are characterized by the thermal mass, width and 
height of the peak of the quasiparticle. In the NLO, the 
thermal-mass correction gives the contribution of order 
\/A, which is nonanalytic in A. [22]. The nonanalytic term 
is obtained from a phase space integral with an infrared 
enhancement of the Bose-Einstein distribution function 
f(E) ~ T/E at small E. Since p(k) appears with f(k) in 
C, S, and the inner product (see, e.g., Eq. (Bll),) con- 
sider the following integral to see the nonanalyticity of 
A: 



d 4 p 

H 
12 



p(p)f(p) 



l + Oi 



In 



rriT \ 



2nT 



)H( 



rriT \ * 



where a± = 6, a 2 = 3 — 67^ + 6 In 2 and a' 2 = 3 
with the Euler constant ^ E = 0.577216. We have as- 
sumed that the spectral function has a form p(k) = 
(27r)sgn(fc°)<5(fc 2 — mf,) to obtain Eq. (88). The nonana- 



lytic term mj- appears in RHS of Eq. (88), although the 
integrant is a function of m 2 ^, which is obtained as [31] 



AT 2 

"24" 



1 - 3 



A 
24tt 2 



(89) 



The thermal masses in the leading (~ A) and the next 
leading (~ A 3 / 2 ) orders are given by the one-loop dia- 
gram and the ring diagram resummation, respectively. 
Therefore the contribution of the phase space integral in 
the NLO gives that of order y/X, and the coefficient is 
estimated as c\ ~ 0.104 [22]. We note that this term also 
contains the next-to-next-leading order (NNLO) correc- 
tions of orders A In 1/ A and A to the shear viscosity from 
Eqs. (88) and (89). One might worry about corrections 
from the thermal width of the spectral function; however, 
this is in even higher-orders because the width is of order 
A 2 T, which is smaller than AT. 

Finally, we identify diagrams contributing to the shear 
viscosity in the NNLO corrections. To our knowledge, 
this is the first identification of the diagrams contribut- 
ing to the NNLO corrections. As noted above, the ther- 
mal mass correction to p{k) is one of them. The correc- 
tions in the NLO to S and C can generally contribute to 
the transport coefficient in the NNLO. However, for the 
shear viscosity, the correction to S does not contribute 
in the NNLO. The reason is as follows: The diagrams 
which contribute to the vertex corrections up to of order 



12 



FIG. 9. Diagrams of order A 3 contributing to T (A) %{k, k'). 




FIG. 10. Collision term at the leading order corresponding to 
3 — > 3 scattering. 

A 2 are shown in Fig. 8. The second and third diagrams 
corresponding to the NLO corrections to S are ring dia- 
grams, which do not depend on the external momentum 
k. However, the contributions of the ring diagrams van- 
ish because the vertex function S for the shear viscosity 
is a second rank tensor, and proportional to k % W — S^k 2 . 
Thus the corrections to S start from of order A 2 that 
are higher orders of the coupling. The diagrams of the 

LO and NLO corrections to z{k,k') are shown in 

Fig. 9. These corrections correspond to quantum and 
thermal loop corrections of order A 3 to the squared scat- 
tering amplitude of 2 — > 2 process, which do not include 
multiple scattering such as 3 — > 3 and 2 — > 4. A typi- 
cal diagram corresponding to 3 —¥ 3 process is shown in 
Fig. 10, which process contributes to the next-to-next-to- 
next leading order, since the squared amplitude is pro- 
portional to A 4 . In summary of the NNLO correction to 
the shear viscosity, the corrections of the thermal mass 
up to of order A 3//2 and the squared scattering amplitude 
of order A 3 contribute to the NNLO, while the corrections 
to S do not. 



V. CONCLUSION AND OUTLOOK 

Wc have given a formulation of a resummation method 
for computing transport coefficients in relativistic quan- 
tum field theory, using the 4 theory as an example. 
We separated the diagrams into the dominant part from 
the other diagrams and reformulated them in a way 
which makes the included physics instructive, by adapt- 
ing Eliashberg's method to a relativistic case in the real- 
time formalism. Our result is summarized to Eqs. (64), 
(73), and (74). The self-consistent equation of the vertex, 
Eq. (73), has a meaning of a kinetic equation, and has a 
form similar to the linearized Boltzmann equation. In the 
leading order of the coupling constant, we recover the ki- 
netic equation of previous works [13] (see Appendix B.) 
The higher-order corrections beyond the leading order 



are systematically incorporated in the formalism for the 
first time and found to be nicely summarized as the 
renormalization of the vertex correction, Eq (64), the 
spectral function, p(k), and the collision term through 

ImT^' 3(k,k') and Imllji(k). In the higher-orders of 
the collision term, both effects of multiple scattering and 
quantum loops are expressed in a power of the coupling 
constant. Their effects also appear in the vertex correc- 
tion. The higher-order correction is important to see the 
convergence of the perturbation theory at finite temper- 
ature. 

We have identified the diagrams up to the next-to- 
next-leading order corrections in the weak coupling ex- 
pansion of 4 theory. The detailed calculation for the 
higher-order corrections will be discussed in our future 
work [32] . We emphasize that the advantage of our dia- 
grammatic method is to enable us identify diagrams con- 
tributing to the higher-order corrections, which is a dif- 
ficult task in kinetic approaches and other diagrammatic 
methods. 

Although the formalism is developed using the </> 4 
theory in the present work, the generalization of it to 
fermionic theories is straight forward: The Bose-Einstein 
distribution function in the self-consistent equation is 
simply replaced by the Fermi-Dirac distribution function. 
The spinor structure is introduced into the vertex func- 
tion and the collision term. Since the Fermi-Dirac dis- 
tribution function hif(E) does not diverge at E = 0, a 
infrared cutoff introduced in Sec. IV A is not necessary. 

Gauge theory is more complicated. Resummation 
of collinear divergences in addition to the ladder dia- 
grams is necessary to obtain the correct result in cer- 
tain orders of the coupling constant, which is called 
Landau- Pomeranchuk-Migdal (LPM) effect [33]. The 
self-consistent equation summing the collinear singular- 
ities over has the form of the Boltzmann equation with 
the collision term in 1 — > 2 and 2 — > 1 process [17, 20, 34]. 

It is also interesting to apply the diagram method to 
critical phenomena, where the hydrodynamic mode and 
the fluctuation of order parameter play an important 
role. We have to take into account these modes in addi- 
tion to quasiparticle modes. 
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Appendix A: Several Identities in the real-time 
formalism 



In this Appendix, we show several useful identities in 
the RJA basis [27]. We start with the standard basis. 
The largest and smallest time equation shows 

Y / r i^-in(PuP2,--- ,Pn) I] e ~ aPl = < ( Al ) 
{ik} U=2 

and Kubo-Martin-Schwinger (KMS) relation shows 

^r ili3 ... in (p llP2l --- , Pn ) J] e (^-K = . (A2) 



{ik} 



i k =2 



By diagrammatic analysis, one finds [27] 



r* li2 ... t (pi,P2,--- ,Pn) 



-I\ li ^M^---,Pn)\{e-^-^. (A3) 



ik=2 



Equations (Al) to (A3) are basic identities in the real- 
time formalism [35]. In the R/A basis, Eq. (Al) becomes 

n 

X] II Uam im(Pm) Y[ e " aPkr aia 2 -a n (Pl,P2, - ■ ■ ,Pn) 
{ik} m=1 ik=2 



n 



m— 1 



x r n 



XPl,P2, ■ ■ ■ ,Pn) 



r R R...R(j>l, P2 , ■ ■ ■ ,p n ) J] ( e ^" 1 - 1) = 0. 

rn—l 



(A4) 



Since e^™ -1^0, 

rRR,..B.(Pl,P2, ■ ■ ■ ,Pn) = 0. (A5) 

Similarly, Eq. (A2) becomes in this basis 

r AA ... A (pi,p 2> - • • ,Pn) = 0. (A6) 

Therefore, the vertex function vanishes, when all indices 
are the same. This means particles are not produced 
from the thermal equilibrium system nor absorbed into 
the system, similar to the vacuum at zero temperature 
and zero density. 

Equation (A3) is reduced to 

r a ia2 --- a JPl:P2, ■ ■ ■ ,Pn) 

n 

= r* li2 ... in (pi,p 2 ,--- ,p n ) n U lm am {p m ) 

m=l 

= -r &1 &2-a n (pi,P2, ■■■ ,Pn) 

x n f(p°m) n 



(A7) 



a m =R 



-f(-P° m ) 



For a two-point function, 

n R (-k,k) = r RA (-k,k), (A8) 

implying D R (k) = D A (k). For a three-point function, 

r RAA(P, 9; k ) = -r A RR.{p, q, k ^ f( p + h) - f(k) ' ^ 
Therefore, 

- (f{P + k) - f(k))r RAA (p,q,k) = r ARR (p,q,k) . 

(A10) 

For a four-point function, 

Tl ARR (pi,P2,P3,Pi) 

r ( , f(p° 3 )M) (ah) 

= -T RRAA (pi,P2,P3,P4)^j — oTTT — o\ ■ 

f(-Pl)f(-P2) 



Appendix B: Relation between Boltzmann equation 
and Kubo formula in the leading order 

Here we explicitly show that Eq. (77) is equivalent to a 
linearized Boltzmann equation in the leading order. This 
was shown in Ref. [13] in a different formalism. Let us 
start with a Boltzmann equation, 



2k^f(x, k) = -C[f] 



(Bl) 



where f{x,k) denotes the distribution function, and C[f] 
is the collision term. In the leading order, the collision 
term has the form of 2 — > 2 collision, 

_ 1 f d 3 k' 1 f d 3 q 1 f d 3 q' 1 



2 J (2tt) 3 2E k , J (2tt) 3 2E q J (2tt) 3 2E q , 
X \M\ 2 (2TT) 4 S^(k + q-k' -q') 
X [fkf g (l + fk>)(l + f g >) 

-(l + /fc)(l + /«)/*'/«'] > 

(B2) 

where M. is the scattering amplitude; Ai = iX for <p A the- 
ory in the leading order. At (local) thermal equilibrium, 
the distribution function has the form 



fo(x,k) 



1 



exp[f3(x)u ti (x)k^) — 1 ' 



(B3) 



where u M (x) is the local velocity, and fc M is the four mo- 
mentum at on-shell, k 2 = to 2 , with mass m. Here we 
consider the shear flow, then the LHS of Eq. (Bl) be- 
comes 

2k»d,f (x,p) = -pf (x,p)(l + / o (z,p))r^*0^ii(z) , 

(B4) 

where J^(fc) = 2k i k j - 2<^fc 2 /3, and 



aij = diUj(x) + djUi(x) - -8 tj d k u k ■ (B5) 
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We have considered the local rest frame, = (1,0). 
Note that at the local rest frame diUj(x) is nonzero, al- 
though u.i(x) — 0. By linearizing the Boltzmann equa- 
tion around the thermal equilibrium f(x, k) = fo(x, k) + 
5f(x, k), one hnds 

C[/]^(l + /o(^))^/0^-(l + /o(^)) 

X / (2703 2E k M Ek,) l Uf^ MEq ' ] 
x (2ir) 4 5^(k + q-k' -q') 



2E k (3f {E k )(l + h{E k ))C^ t ^-a l3 , 



(B6) 



where fo(E q ) = fo(x, k) at the local rest frame, ^Boitz is 
the linearized collision operator, and we chose 8f(x, k) = 
Pfo(E k )(l + f (E k ))^(k)a lJ /2. Then, the linearized 
equation reads 



r^{k) 
2E k 



This is the same form as Eq. (76). Let us check 
C = CboUz. at weak coupling. For this purpose, let us 
estimate the four-point function, ImT*-" 4 ' s(k,k'), and 
Im7T] ? wo " loop (fc). Using Feynman rule, Eqs. (28) to (36), 
we find 



ImT (A) 3 (M') 
1 

~ ~2 



-f{q'°)p{q'){2^8{k + q-k' -q 1 ), 



d 4 q' 
(M 



The collision term becomes 
C^{k) = -2 -Ll m 77^- 1 °°P(fc)^(fc) 



i r d 4 k' 



2k° J (2tt) 4 

f-Hk^x 2 r d 4 k' 



2fc° 6 J (2tt) 4 



P {k')lviT {A) z (k,k')^{k') 
f(k'°) P (k r ) 



(l + /(9°))p(?) / ^fW°) P W) 



(2tt) 4V j ^ J (2^)4- 

x {2n) 4 5 (4) {k + q-k' -q') 

x (k)-^(k')+ <p ij {q)-<p iS (<?')]• 

(Bll) 

Here we have used the relation ip l ^(—k) = —^{k) de- 
rived from Eq. (75) to obtain the third line. Equa- 
tion (Bll) includes positive and negative energies, while 
Eq. (B6) includes the only positive energy. This colli- 
sion term contains three processes: 1 — > 3, 2 — > 2, and 
34 1 scatterings. 1 — > 3 and 3 — > 1 collisions can be 
neglected in the leading order because the on-shell condi- 
tion of quasiparticle does not satisfy. Then, the collision 
term for the positive energy state, k° > 0, becomes 



( B7 ) &p ij (k) 



f Hk ) A 2 f d*k' /(fc /o )p(A ./ )e(fc / 0) 



2k° 2 J (2tt) 4 
d 4 q 



(2tt) 4 

d 4 q' 



(l + f(q°))p(q)6(q ) 
./V>(<z')%'°) 



(27T)*- 

x {2^) 4 5 {4 \k + q-k' -q 1 ) 

x [^(k)-^(k') + ^(q)-ip»(q')]. 

(B12) 

By using the quasiparticle approximation, 

P(k)9(k°) = ^(27T)S(k° - E k ), (B13) 
(B8) with E k = yjk 2 + m^, we obtain the collision operator 



and 



ImTI 



two-loop 
R 



(k) 



A 2 /■ d 4 q 

12 y (27rp' 



^ fe °-i)^ / ^(l + /(g°))p(?) 



x (2^) 4 <5 (4) (fc + g-fc'-g')- 



Comparing Eqs. (B8) and (B9), we find 
1 f d 4 k! 



ImH 



two-loop 

R 



(k) 



(B9) 



6J (2vr) 



p{k')lmT {A) 3 (k,k') 



(BIO) 



£^'(fc) 

f^ 1 (E k ) A 2 



d 3 fc' 



2£ fc 2 7 (2tt) 3 2£* 
d 3 9 1 



/(^) 



(2tt) 3 2£ 9 v J v 9 ^ y (2tt) 3 2£7,» 

x (27r) 4 (5 (4) (fc + q-fc'-q') 

x [<p ij (k) - <^'(fc') + ip^{q) - ip ij (q')] . 

(B14) 

Thus, the collision operator in the Boltzmann equation is 
equivalent to that in diagrammatic method in the lead- 
ing order, L = ^Boitz.- In the neutral scalar theory, the 
negative energy state is identical to the positive energy 
state, so that the collision term for the negative energy 
state is the same as that for the positive energy state. 
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Appendix C: Some properties of the collision 
operator 

The collision operator of the linearized Boltzmann 
equation has some basic properties: semi-positive defi- 
niteness, self-adjointness, and the conserved charges be- 
ing the collision invariants. Here, we focus on the self- 
adjointness and the conserved charges being the collision 
invariants, although the semi-positive definiteness of the 
collision operator is necessary to ensure the stability of 
the thermal equilibrium. 



1. Self-adjointness of the collision operator 

In order to show the self-adjointness, we write the lin- 
earized collision operator as 

cm = / 0f w{k ' )c{k > k>{k,) ' (C1) 



where 



C(k, k') = — lmn R (k) W^tk'Wk - k') 



+ ^lmT^ 3 3(k,k')W-\k')p(k') . 

(C2) 

This is the momentum representation of the linearized 
collision operator. The Hermite conjugate is defined as 



which means that 

<■ d 4 k 



(C3) 



(2tt)4 

d 4 k' 



(27T) 



W{k)y\{k) 
^W(k')C(k,k')Mk') 



d 4 k 
d A k' 



W(k)C*{k,k')<px{k) 



(2tt)4 
d i k 



W(fc') 

W(k)C^(k',k)ipi(k) 



or 



tf(k',k)=L*(k,k') 



(C4) 



(C5) 
(C6) 



Let us check the self-adjointness or Hermiticity, & = C. 
The first term of RHS in Eq. (C2) is obviously symmetric 



under k O k' . The four-point function, 3 (fc,fc'), in 

Eq. (C2) has a symmetry, 

from Eq. (All). Thus the second term of RHS in 
Eq. (C2) is symmetric: 

imr^ 3 3(fc,fc')r 1 (fc'°)(i + /(fc' ))~ 1 
= imr(- 4)3 3(fc',fc)r 1 (fc°)(i + /(A : ))- 1 . 

Therefore, the collision operator is Hermite, = C. 
More precisely, the collision operator is symmetric be- 
cause Eq. (C2) is a real function. 



2. Collision invariant and Ward-Takahashi identity 

In this subsection, we derive Ward-Takahashi (WT) 
identity in the R/A basis, and show the WT identity 
implies that the conserved charges are zero modes of 
the collision operator, Eq. (77). In the imaginary-time 
formalism, WT identity is derived in Refs. [17, 36, 37]. 
The WT identity is useful to constrain diagrams in or- 
der to keep symmetries in the resummed perturbation 
theory [17]. Conversely, to obtain the correct result in 
a certain order of perturbation theory, the collision term 
must satisfy the WT identity. 

The conserved current in <f> A theory is a energy- 
momentum current, defined by 



(C9) 



In particular, the charge density of the momentum cur- 
rent is the bilinear operator, 



Ji 



di<f>d°<P . 



(CIO) 



This bilinear form is important to relate the WT identity 
and collision term, as we will see later. In order to derive 
the WT identity in thermal bath, consider the following 
correlation function: 



(T c j?(z)<p(x)ct>(y)} , 



(Gil) 



where Tc denotes a complex-time-path ordering shown 
in Fig. 1. Taking <9 2M , we obtain 



= d c (z°-x o )(r c [j°(z),0(x)]^)) 



(C12) 



where we used the conservation law dpj? = 0. Noting 
that the equal-time commutation relation gives an in- 
finitesimal translation. 



\j° i {z),4>{x)] = -i8 ( ?\z-x)d i( l>{x) 



(C13) 
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we obtain 

d ztt {Tcj?(z)<Kx)<l>{y)) = -i8^\z-x){Tcd i( l>{x)<j>{y)) 

-iS^(z-y)(T c( p( X )d i( t>(y)). 

(C14) 

This is the WT identity in the coordinate space on the 
complex time path. In the real-time formalism, the com- 
plex path is decomposed into four parts as shown in 
Fig. 1. In the standard basis, Eq. (C14) becomes 

d^{T c % c {z)r{x)<t> b {y)) 

= -iX cd a 5^(z - x){Td zi cj> d (z)<t> b (y)) (C15) 
-iX b d 8^{z-y){T<j> a {x)d zi <p d {z)), 

where j° c (x) = X cde di^ d (x)d°(jf{x) with X cde = 
$cd9de- Using the propagator and the three-point vertex 
r M (z, w, w'), we obtain 



d^j d A w J d 4 w'r^ de ( Zl w,w')D da (w,x)D eb (w\y) 

= -iX cd a S^(z-x)d zi D db (z,y) 
-iX b d 8^\z-y)d zi D da {z,x). 



In momentum space, this becomes 
ip^ cde {p,q,k)D da (q)D eb (k) 



hX cd a D db (k)-q t X c h d D da (q), 



(C16) 



(C17) 



where p, k and q are not independent but satisfies 
p + k + q = 0. Multiplying both sides of Eq. (C17) by 
D- a \(q)D^(k),wc£nd 

W^ ca , b ,(p 7 q, k) = -hX cb , d D-^{q) - qi X c d a ,D- b )(k) . 

(C18) 

At tree level, this is reduced to 



iP^mcab&V^) = ~ k i X cbD- F \ da {q)-q % X c d a D- F \ cb {k), 

(C19) 

where F^^.^ij), q, k) is the tree-level vertex, 

r (0)i;cab(P> <i> k ) = Xcab(-(q t k» + k^) + <5f (k ■ q + TO 2 )) , 

(C20) 

and D F b (k) is the free propagator. Here, we decompose 
the full vertex-function into that at tree level and its 
correction, Sr^ abc : 

r» abc ( Pl q , k) = r[ mabc ( P , q , k) + sr» abc ( P , q , k) . (C2i) 

Then, the correction to the vertex satisfies 

iPij.dr^ cab (p,q, k) = -ikiX cb d II da (q) - iq.X^ a II db {k) . 

(C22) 

In the vacuum, Eq (C22) gives the relation between the 
charge and wave function renormalization factors. In ad- 
dition to that, at finite temperature, we will show that 



conserved charges are the collision invariants of the colli- 
sion operator. To see this, we, first, take the limit p — >• 0, 

woSr° cab (p, q ,k) = -ik t (x cb d n da ( q ) -x d a n db {k)) 

= F. llcab {k,p°). 

(C23) 

As we mentioned above, jf is a bilinear operator, so that 
the three-point vertex can be written by using the four- 
point vertex T a bcd as 

f d 4 k' 

sr °cab(P,q, k ) = j j^ri(-i)T a bde{q,k,-k',-q') 

xDf d (k')D^q')r^ i . cgf (p,q',k'), 

(C24) 



so that Fi. cab (k,p°) becomes 



Fi-cab(k,p°) = -i 



d A k' 



T a bde(q,k,-k',-q') 



(27T)* 

xk , i (X c d g D°°tf)-X ef e Dt d (kr) 

-Df d (k')D^(k')F i;cgf (k\p°)). 

(C25) 



Let us define T (A) a fcde(<7, k, -k', -q') by 

Tabde(q, k, -k' , -q') 

= T (A) ab de(q, k, -k', -q') 
d 4 k" 



j^ i {--'i)Tabd'e'{q,k,-k ,-q ) 

xD b ' d '(k")D a ' e \q")T (A) a ^ de (q",k",~k',-q'). 

(C26) 

This is equivalent to Eq. (57) but in the standard basis. 
Comparing Eq. (C26) with Eq. (C25), we find 



Fi- lCa b(k,p°) = - i 



^^Ws^-fcW) 

(A) 

x (X d g D^{q') - X cf e D^(k')) . 

(C27) 



We have derived Eq. (C27) in the standard basis; how- 
ever, this is independent of a choice of basis. Now, let us 
apply it to the R/A basis. We take the case a = b = A, 
and c = R: 

X RA R (p,q,k) = l, (C28) 
X RA A (p,q,k) = l + f(p°) + f(k ), (C29) 
X RR R (p, q,k) = l + f(p°) + f(q°) , (C30) 
X RR A {p,q,k) = 0. (C31) 
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Therefore, F^raa becomes 

Fi.RAA(k,P°) 

\t^aarr(-P - k, k, -k',p+ k') 
(2tt) 4 L 

x {-i)(D R {p + k')-D A (k')) 
+ T {a) aara(-P - k, k, -k',p + k')f(k°) 

x (-i)(D A (p + k')-D A (k')) 
- T {a) aaar(~P - k, k, -k',p + k')j{k Q + p°) 

x (-i)( D R (p + k')-D R (k'))] , 

(C32) 

where we have omitted 1 + f(p°) in the second (third) 
term because there is no pole in the lower (upper) com- 
plex k'° plane. Taking the limit p a — > 0, we find 

lim Fi.RAA(k,p°) 

Po->0 



d A k' 



KT {A) AAR R (~k,k,~k' \k') P {k') . 

(C33) 



On the other side, from Eq. (C23), 

Fi- R AA(k,p°) = -ik l (X RA R {p,q,k)IlR A (q) 

-XR R A (p 7 q,k)IlR A (k)) (C34) 

= -iki(n A {q) - n A (k)) . 



Taking p° ->• 0, we find 



lim Fi. R AA(k,p 



= 21m n R (k) h 



(C35) 



Inserting Eq. (C35) into Eq. (C33), we find 



-2Im77 fl (fc)fc J + J ^ImT^ 3 3(fc,fc>(^ = 0, 

(C36) 

where Im T (A)3 3 (fc, k') = lmT {A) AARR (-k, k, -k' , k') is 
used. Using the collision operator, Eq. (77), we find 



Ch = . 



(C37) 



Therefore, the conserved charge ki is zero mode of the 
collision operator. 
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